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We report the observation of a large linear magnetoresistance (MR) and Shubnikov-de Hass (SdH) quantum 
oscillations in single crystals of YPdBi Heusler topological insulators. Owning to the successfully obtained the 
high-quality YPdBi single crystals, large non- saturating linear MR of as high as 350% at 5K and over 120% at 
300 K under a moderate magnetic field of 7 T is observed. In addition to the large, field-linear MR, the samples 
exhibit pronounced SdH quantum oscillations at low temperature. Analysis of the SdH data manifests that the 
high-mobility bulk electron carriers dominate the magnetotransport and are responsible for the observed large 
linear MR in YPdBi crystals. These findings imply that the Heusler-based topological insulators have 
superiorities for investigating the novel quantum transport properties and developing the potential applications. 

Materials exhibiting large linear magnetoresistance (MR) have attracted intense research interest due to 
their potential applications in magnetic random access memory and magnetic sensors \ The MR of non- 
magnetic metals with open Fermi surfaces {e.g., Au) can be linear and non- saturating at high fields^. 
Large linear MR effects have recently been explored in materials with approximately zero bandgaps, such as 
doped silver chalcogenides^ single -crystal bismuth thin films^, indium antimonide^, and multi-layer graphene^'^. 
There are two prevailing models that explain the origin of linear MR behaviors. One is the quantum model by 
Abrikosov^ for materials with gapless linear dispersion spectra and the other is the classical model by Parish and 
Littlewood (PL model)^ for strongly inhomogeneous systems. The recent discoveries of giant linear MR in 
Bi2Se3^°'^^ and Bi2Te3^^'^^ topological insulators (TIs) has brought renewed interest to linear MR. Because of their 
unusual surface states that are naturally zero bandgap with linear dispersion, TIs provide a perfect platform on 
which to study the origin of linear MR. Because some Heusler compounds have recently been predicted to be 
TIs^^"^^, it is worth investigating whether large linear MR can also be observed in Heusler alloys. Considerable 
efforts have been devoted to the study of the electric structures and transport properties of a number of Heusler- 
based TIs^^"^^. X-ray photoelectron spectroscopy results have shown that such Heusler-based materials as YPtSb, 
LaPtBi, LuPtSb and LuPdBi are gapless semiconductors with very high mobility^ ^'^^'^^ although large linear MR 
and desired Shubnikov-de Haas (SdH) quantum oscillations have not been observed experimentally in these 
materials. Among these predicted Heusler TIs is YPdBi, which stands the border between the trivial and topo- 
logical insulator states. Here, we present experimental evidence demonstrating the existence of large linear MR 
and the first observation of SdH quantum oscillations in high-quality single crystals of YPdBi. The MR is about 
350% and 120% at 5 and 300 K, respectively, under a magnetic field of 7 T. Analysis of temperature- and angle- 
dependence of SdH oscillations manifest that the observed large linear MR features do not come from the surface 
states, but arise from the high-mobility 3D bulk electron carriers in YPdBi crystals. 

Results 

Structural characterizations. High-quality single crystals of YPdBi were grown by the flux method as described 
in the Methods section. Typical dimensions of the single crystals were about 1.0 mm X 1.0 mm X 1.0 mm and 
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Figure 1 | (a) Photograph of a single crystal of YPdBi. (b) Laue diffraction pattern of sing-crystal YPdBi generated with the beam axis coniccident with the 
[100] zone axis at room temperature, (c) Powder XRD pattern of pulverized single- crystal YPdBi sample, together with the simulated XRD data for 
fully ordered YPdBi with MgAgAs-type half-Heusler structure. Left inset shows a structure view of half-Heusler and right inset shows the X-ray rocking 
curve of the as-grown YPdBi single crystals, (200)- Bragg reflection. 



the dimensions of the largest crystals were approximately 5 mm X 
4 mm X 4 mm as shown in Figure la. Most of the crystals were 
located at the bottom of the crucible, which could have been 
caused by the high density of the materials. The quality of the 
single crystals was examined by the x-ray Laue back- reflection 
technique from different positions across one facet. The resulting 
Laue image obtained from a piece YPdBi single crystal, shown in 
Figure lb, clearly indicates that the sample was indeed a single crystal 
with very well- developed (100) facets. Small triangular-hke steps 
were found on the surfaces, which strongly suggests that the 
crystals grew through two-dimensional nucleation and a layer-by- 
layer growth mechanism on the dense planes. The composition of 
the crystals was charcaterized by energy- dispersive x-ray analysis 
(EDX) (see supplementary material Figure SI and Table SI). The 
EDX results revealed that the average ratio of Y to Pd to Bi in the 
crystals was about 34.2 to 33.3 to 32.6, which was very close to to the 
stoichiometric ratio. The uncertainty of measurement was less 1.0%, 
which was determined by repeating the measurements on the same 
sample surface. No flux component was detected. To deter- 
mine the crystal structure of the samples, we ground the crystal 
and performed power X-ray diffraction (PXRD) measurements. 
Based on the PXRD patterns shown in Figure Ic, it is evident that 
these crytsals were single-phase crystals. The structure was found to 
be a MgAgAs-type cubic crystal structure (inset of Figure Ic) with 
lattice parameters of 6.68A, which was consistent with previous 
reports on polyscrystalUne samples^^"^\ Moreover, the normal axis 
of the crystal platelets was along the (100) direction as determined by 
Laue diffraction. The quality of the crystal was further examined 
using the X-ray rocking curve, as shown in the inset of Figure Ic. 
The typical full-width at half-maximum (FWHM) of the X-ray 
rocking curve for the (200) reflection was about 0.12°, which again 
indicated that the crystals had good crystallinity and that the growth 
technique was suitable for this type of alloy. 
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Temperature-dependent resistivity. In Figure 2a and b, we show the Figure 2 | (a) Temperature-dependent zero-field resistivity p^x- In the 

temperature dependence of the zero-field resistivity, pxx, and the region of constant carrier concentration, 5 K< r< 100 K, the resistivity 

carrier concentration derived from measurements of the Hall decrease with decreasing temperature owing to decreased phonon 

effect, assuming a one-band model. The zero-field resistivity curve scattering, (b) Variation of carrier concentration with temperature 

shows a broad maximum at T ~ 100 K and a weak temperature deduced from Hall effect measurement. The inset shows the Hall resistivity 

dependence over the whole temperature range, i.e., the dp/dT was pxy of YPdBi as a function of magnetic field. The dominant charge carriers 

quite small. Below the peak, the behavior of the resistivity was are electrons. 
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metallic-like, where the electron carrier density, n^, was mostly 
temperature independent. In this constant carrier density region, 
the resistivity decreased with decreasing temperature owing to 
decreased phonon scattering. Above the peak, the resistivity 
decreased and the carrier density, ng, increased with increasing 
temperature, which is commonly found in semimetal or gapless 
semiconductors that do not have large activation energies. The 
magnitude of pxx (~0.8 mQ cm at 300 K) is consistent with the 
reported value of the polycrystalline YPdBi compound^^. The high 
value of the carrier density {n^ ~ 2.5 X 10^^ cm"^ at 300 K) indicated 
that YPdBi is a high-carrier gapless semiconductor. 

Figure 3a shows an overview of the temperature dependence of the 
resistivity, p(T), of YPdBi measured under a series of appUed mag- 
netic fields. The field was applied perpendicularly to the current 
direction. Strikingly, the resistivity increased dramatically over a 
wide range of temperature, from 5 to 300 K, when a 5 T magnetic 
field was applied. Another interesting feature of the 5 T curve is that 
the overall behavior of the curve changed, e.g., the broad peak 
disappeared, a plateau formed below 100 K and the temperature 
dependence became much stronger, i.e., the dp/dT was much larger 
than in the zero -field curve. In other words, the ratio of the peak value 
to the value at 300 K, p(100 K)/p(300 K), was about 1.3, whereas it 
increased to 2.1 in the 5 T curve. All these changes under a moderate 
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Figure 3 | (a) Temperature- dependent resistivity of the single -crystal of 
YPdBi at a series of magnetic fields H = 0, 1, 2, 3, 4 and 5 T applied fields, 
(b) Normalized MR as a function of magnetic field H at a series of 
temperatures T = 5, 50, 100, 200 and 300 K. Here MR = [p(H) - p(0)]/ 
p(0)] X 100%, where p(H) and p(0) are the resistivity with and without the 
magnetic field H, respectively. Inset, the observed MR of YPdBi at 300 K 
(open circles) and fit (black line) with the quadratic equation MR = a\H\ 
+ bH^. 



magnetic field (5 T) indicate that the electron behaviors were dra- 
matically altered by the magnetic field. 

Linear magnetoresistance. Figure 3b shows the magnetic-field 
dependence of the normalized MR at different temperatures. The 
MR is defined as [p(H) - p(0)]/p(0)] X 100%, where p(H) and 
p(0) are the resistivities at field H and zero field, respectively. A 
robust MR, as large as 350%, was observed under a 7 T magnetic 
field over a broad temperature range from 5 to 100 K. More 
importantly, the MR was temperature independent over this tempe- 
rature range and near linearly proportional to the magnetic field. 
These two characteristics are essential for application of high 
magnetic field sensors at low temperatures in a wide temperature 
range. The field sensitivity was about 50%/T. Above 100 K, the MR 
gradually reduced with temperature, and reached to 120% at 300 K, 
which is also remarkable for field sensors at room temperature. 

To understand the mechanism/or physics behind this giant MR, 
we also measured the MR when the magnetic field was parallel to the 
current direction, or longitudinal MR (see supplementary material 
Figure S2). We found that although the longitudinal MR showed 
non- saturating behavior, its magnitude was relatively smaller than 
that of the transverse MR. At temperatures below 100 K, the longit- 
udinal MR had a weak quasi-linear dependence on the applied mag- 
netic field and it was almost independent of temperature. At 
temperatures higher than 100 K, however, the quadratic contri- 
bution in the low fields became more pronounced. The longitudinal 
MR of the polycrystalline sample was also non -saturating, but only a 
very small MR of 2% was observed under a magnetic field of 7 T (see 
supplementary material Figure S3). 

The inset of Figure 3b shows the best fitting result for the MR curve 
at 300 K. As a result, we found that the relation between MR and the 
magnetic field, H, could be described by the following quadratic 
equation: MR = a\H\ + bH^. From the above observations, we can 
conclude that the MR originated from contributions of both the 
linear and parabolic terms. The parabolic term is well understood 
as a result of the Lorentz force, whereas the origin of the linear MR 
remains controversial. To ascertain the origin of the observed linear 
MR, we first considered the quantum linear MR description pro- 
posed by Abrikosov^ in which the zero-gap band structure with linear 
dispersion is important. One thus would expect that the linear MR 
observed in the YPdBi single crystals would be explained by the 
Abrikosov's quantum MR^, as the bands of Heusler compounds 
exhibit linear dispersion close to the Fermi energy and the charge 
carriers behave like massless Dirac fermions. Therefore, the MR 
should be linear from low to high field and it should be independent 
of temperature. It is almost independent of temperature for T < 
100 K for magnetic fields higher than 2 T. In addition, if the linear 
MR was quantum mechanics in nature, then as the temperature 
increased, we would expect the MR to increase because the linear 
energy dispersion of the topological surface state can persist at both 
low and high temperatures. The linear MR in epitaxial grapheme^ 
and topological insulator Bi2Te3 nanosheets^^ was found to increase 
with increasing temperature and was therefore interpreted as due to 
Abrikosove's linear quantum MR^. However, there is clear evidence 
that the linear MR in our crystals decreased with increasing temper- 
ature for T > 100 K. 

Shubnikov-de Hass oscillations. Shubnikov-de Haas (SdH) quantum 
oscillations has been identified as a convincing tool for characterizing 
quantum transport in materials showing the three-dimensional (3D) 
bulk^^'^^ and the two-dimensional (2D) surface states^^"^^, respec- 
tively. We therefore carried out low-temperature and high-field 
magnetoransport measurement to provide experimental evident for 
the quantum transport in YPdBi crystal. The magnetic field is 
perpendicular to both the current flow and the surface of the 
YPdBi crystal. The magnetic-field dependent resistivity R^ shows 
traces of SdH oscillations in our raw data (See supplementary 
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Figure S4), reflecting the high carrier mobility. In Figure 4a, we show 
the oscillatory component of AR^ versus 1/H at various tempe- 
ratures (T) after subtracting their linear background (See supplemen- 
tary Figure S4). The amplitude of the SdH oscillations decreases with 
increasing T, and the oscillations are not observed for T > 20 K. 
Moreover, a single oscillation frequency can be extracted from fast 
Fourier transform (FFT) spectra (/sdH(T) = 45 T, See supplementary 
Figure S5). The frequency directly gives a Fermi surface cross-section 
area of Sp = 0.037A"^. The tilt angle between H and the surface of 
crystal can be varied from 90° to 0°, and we find the SdH oscillations 
can be also observable at 0° (See supplementary Figure S6), sug- 
gesting the oscillations are essentially due to the Landau quanti- 
zation of the 3D bulk Fermi surface. Assuming a spherical Fermi 
surface, the carrier concentration at 2 K was estimated to be 1.7 X 
10^^ cm"^, which is close to that estimated by the Hall analyses, so we 
can ascribe the origin of the SdH oscillation not to the possible surface 
states but to the intrinsic gapless bulk nature of YPdBi. In addition, 
the temperature dependence of the SdH oscillation amplitude can be 
fitted to the standard Lifshitz-Kosevich theory^ \ 

ARiT,H)^ eM-2n'KnTo/AE^{B)] . ltY^rf^^\^^ 

In Eq. (1), A£^ and are the fitting parameters, and H is the 
magnetic field position of the Nth minimum in R. IS.E^{B) = 
heH/lnnf is the energy gap between the Nth and {N + l)th LL, 
where is the effective mass of the carriers, n is the electron charge, 
and h is the Planck constant. TD = h/ Au^tkb is the Dingle 
temperature, where t is the quantum lifetime of the carriers due to 
scattering, and Kb is Boltzmann's constant. In Figure 4b, we show the 
temperature dependence of SdH oscillation amplitude at 9 T. The 
solid line shows the best fit to Eq. (1). The inset of Fig. 4b shows 
the calculated A£jv as a function of H. The slope of a linear fitting 
yields a rather small cyclotron mass ~ 0.09 (m^ = 9.1 X 
10"^^ kg is the mass of free electron). The Dingle temperature is 
found to be 6 K from the slope in the semilog plot of 
^R/R{0)H sm\i{2Ti^ KbT / ^E^) vs. 1/H at T = 2 K [Figure 4c]. 
From To we extract the carrier scattering time t = 4 X 10~^^s. 
This value of t corresponds to the carrier mobility of 7800 cm^ 
V~^s~^ at 2 K, in reasonable agreement with the value of the Hall 
mobility, which provides a strong, quantitative argument that the SdH 
oscillations arise from the 3D bulk states. 

Discussion 

To understand the origin of the observed large linear MR, we recalled 
the classical PL modeF, in which the linear MR is expected to be 
governed by carrier mobility. The classical PL model predicts that the 
crossover field. He (the field at which the MR curve becomes linear), 
is in inverse proportion to the carrier mobility, and that it should 
continually increase with increasing temperature due to decreasing 
mobility. In Figure 5a, we show the He and carrier mobility, p.^, as a 
function of temperature. The He was obtained from the first- order 
derivative of the MR curve with field as plotted in the inset of Fig. 3a. 
The carrier mobility, ii^{T) = Rh{T)/p{T), was extracted from the 
Hall coefficient, RhCOj using a single band model. Here, we see that 
He remains somewhat constant for T < 100 K, whereas it continu- 
ally increases with increasing temperatures for T > 100 K. Moreover, 
we found that the Hall mobility, fi^, decreased slightly with temper- 
ature, although the magnitude was not strongly dependent on the 
temperature. The PL modeP also predicts that the slope of the linear 
part of the MR, i.e., dMR(T)/dH, is proportional to the carrier mobil- 
ity, ^e- Figure 5b displays a plot of dMR(T)/dH as a function of ^e- 
Interestingly, we found that dMR(T)/dH depends crucially on carrier 
mobility. The higher the mobility, the larger the dMR{T)/dH mag- 
nitude. We suggest that MR and fi^ qualitatively follow the relation 
MR(T) = iiXt)H, which has been reported in some recent works on 
polycrystalline samples^""^^. We can therefore conclude that the 
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Figure 4 | SdH oscillations in YPdBi crystal, (a) The variations of AR^ vs. 
\IH, after subtracting a fitted smooth background at various temperatures, 
(b) Temperature dependence of the SdH oscillation amplitude AR^ (T)/ 
R^ at 9 T. The solid line is a fit to the Lifshitz-Kosevich formula, from 
which we extract Landau level (LL) energy gap AE^v- The inset shows A£jv 
as a function of H (magnetic field positions of minima in R^ 
corresponding to different LLs), and the effective mass: —0.09 mg is 
extracted from the slope of the linear fitting, (c) Dingle plots of In 
[Ai^/i^(0)H sinh {2ti^ KbT / ^E^)\ versus 1/H at different temperatures. The 
Dingle temperature T^ = 6 K is calculated from the slope of the linear fit, 
corresponding to a carrier lifetime x~4 X 10"^^ s and an effective mobility 
of 7800 cm^ V"V\ 

linear part of MR of Heusler compounds is mainly controlled by 
mobility, particularly in high temperature regions. According to 
the classical PL MR modeF, linear MR is a consequence of mobility 
fluctuations. Local probe techniques showed that electron-hole 
puddles^^ and charge imhomogeneity^^ were found in high-quality 
graphene in which the Dirac electronic dispersion spectrum is firmly 
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Figure 5 | (a) The crossover field Hq on left axis and Hall mobility fi^ on 
right axis as a function of temperature. Error bars are standard deviation. 
Inset shows the first derivative of MR of YPdBi with field at T = 200 K, the 
intersection of linear indicates a turning point between parabolic and 
linear MR. (b) Linear part of the slope of the MR, i.e., dMRiT)/dH as a 
function of carrier mobility fi^. Error bars are standard deviation. 



established. We therefore cannot rule out the possibility that strong 
electronic inhomogeneity is responsible for the linear MR in our 
YPdBi single crystals. 

Let us summarize the important features of high-quality YPdBi 
crystal. First, it shows a large non- saturating MR as high as 350% at 
5 K and over 120% at 300 K under a magnetic field of 7 T. The large 
linear MR exhibited an interesting dependence on the temperature, 
and the magnitude of the MR increases markedly as mobility 
increases at T < 100 K. Second, it has a gapless semiconductor band 
structure with high-mobility 3D bulk carriers. The high-mobility and 
relatively low carrier concentration are the key factors for the obser- 
vation of SdH oscillations at low temperature. Third, the high- 
mobility electron conduction is rather sensitive to the randomly 
oriented domains, as evidenced by the absence of SdH oscillations 
in the polycrystalline Heusler TIs with ultrahigh mobility carriers^^. 
To reveal the origin of the high electron mobility, we have performed 
band- structure calculations for YPdBi as shown in supplementary 
Figure S7. The results are consistent with previous calculations^^, and 
at the equilibrium lattice constant (a = 6.68A) a linear and gapless 
(Dirac-cone-like) energy dispersion is formed by the Fg hole-like and 
the highly dispersive electron-like bands. Moreover, as pointed 
out by Chadov et a/.^^, a very small variation of the lattice constant 
can further shift down the parabolic dispersion heavy-hole Fs band, 
leaving the single Driac cone at the Fermi energy, which allows all 
conduction carriers to condensate at the lowest Landau level. These 
features are compatible with the observation of the SdH oscillations, 
suggesting an isotropic Fermi surface in YPdBi. Nevertheless, a 
precise direct observation of Fermi surface by angle- resolved 



photoemission spectroscopy (ARPES) is desired to indentify the 
band structures at Fermi surface in YPdBi. 

Methods 

Sample preparation. Single crystals of YPdBi were grown in a two-step process. The 
first step was to prepare polycrystalline samples of YPdBi by arc melting 
stoichiometric amounts of the constituent elements in a high-purity argon 
atmosphere. Y (99.95%) pieces, Pd (99.99%) grains, and Bi (99.999%) grains were 
used as starting materials. An excess of 3% Bi was used to compensate for the loss of 
elements during the arc melting. For better homogeneity and crystallinity, the arc- 
melted ingots were wrapped in Tantalum foil and then annealed in evacuated quartz 
tubes at 1073 K for two weeks. The second step was to grow crystals of YPdBi by a flux 
method. The fabricated polycrystalline YPdBi was grounded, mixed with fluxed Pb 
pieces in an atomic ratio of 1:10, and placed in a Tantalum crucible, which was loaded 
into a fused quartz tube. For estimation of the best ratio of the fluxed material to 
YPdBi, we examined the mixtures with various weight ratios, ranging from 20: 1 to 1 : 1 . 
The 10:1 mixture resulted in the best solubility and crystal growth. The tube was 
sealed under Ar gas at the pressure of lO"'* Torr and then placed in a furnace. We also 
tried to estimate the best cooling rate. The best growth conditions at the 10:1 ratio 
were as follows: the mixed sample was heated at 900 °C and then at 1300°C for 15 
hours each and then slowly cooled to 1 030° C at a rate of 3°C/h. 

Materials characterization. The composition of the single crystal samples was 
determined by energy-dispersive X-ray (EDX) spectroscopy and the structure were 
checked by X-ray diffraction (XRD) with Cu-Ka radiation. The single-crystal 
orientation was checked by a standard Laue diffraction technique. 

Transport measurements. We milled the single-crystal samples to a thickness of 0.20 
±0.01 mm for the transport measurements. Electrical leads were attached to the 
samples using room-temperature cured silver paste by gold wires. The in-plane 
resistivity, pxx, and Hall resistivity of a selected crystal were all measured using the 
four-terminal method with a dc-gauge current of 5000 |iA between 2 and 300 K, using 
a commercial Physical Properties Measurement System (PPMS) from Quantum 
Design. The perpendicular and longitudinal magnetoresistances were measured with 
the magnetic field perpendicular or parallel to the electrical current direction. 
Magnetoresistance versus magnetic field was also measured at several temperatures. 
The Hall effect was measured by rotating the crystal by 180° in a magnetic field of 5 T 
in PPMS between 2 and 300 K. The Hall coefficient was calculated from the slope of 
the measured Hall effect curves. 

Band structure calculations. The electronic structure calculations in this work were 
performed using full-potential linearized augmented plane- wave method, as 
implemented in the package WIEN2K^^. The exchange correlation of electrons was 
treated within the local spin density approximation (LSD A) including Spin-orbital 
coupling (SOC). Meanwhile, al7X 17X 17 /c-point grid was used in the calculations, 
equivalent to 5000 k points in the first Brillouin zone. Moreover, the muffin-tin radii 
of the atoms are 2.5 a.u, which are generated by the system automatically. The lattice 
parameters and atomic positions were taken from our experimental data. 
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